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Abstract 



We illustrate how to extend the concept of structural stability through ap- 
plying it to the front propagation speed selection problem. This considera- 
tion leads us to a renormalization group study of the problem. The study 
illustrates two very general conclusions: (1) singular perturbations in applied 
mathematics are best understood as renormalized perturbation methods, and 
(2) amplitude equations are renormalization group equations. 
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I Introduction 



When a very thin film made of diblock copolymers [0, ^ in the disordered 
phase is quenched sufficiently, microphase separation occurs, and segregation 
patterns are formed. What happens if we cool the film from one end? We 
would expect the appearance of a segregation pattern invading the featureless 
disordered phase. The quenched film in the disordered state is thermodynam- 
ically unstable. Thus to facilitate the observability of such propagating front 
phenomena, the growth of spontaneous fluctuations before the front must be 
suppressed. This could be accomplished, for example, by sliding a cooling 
block along the film. If we slide the block too quickly or too slowly, however, 
we would not observe any intrinsic front invasion behavior into the disordered 
phase; if it is too fast, the unstable phase may spontaneously order before the 
front invasion, and if it is too slow, the invasion is restricted by the presence 
of the cooling front. What is the natural speed, given the quench depth? How 
does the pattern invade with this 'natural speed' ? For example, suppose the 
equilibrium pattern for the low temperature state is a triangular lattice. When 
this phase invades into the quenched disordered phase, do we observe the tri- 
angular lattice immediately, or do we observe a lamellar phase first, which 
later orders into the triangular lattice? What are their speeds? 

Now, let us examine an example. Perhaps the simplest model of diblock 
copolymer melt dynamics is the following partial differential equation [^, ^, |^ 

dtip = A{-Ti; + gip^ - DAip) - B{i) - a), (1.1) 
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where ijj is the order parameter field, r, g, D and B are positive constants, 
and a is a constant which could be negative. Figure 1 illustrates the quench- 
ing process due to the moving cooling front simulated by the cell-dynamical 
system]^, |Tl| corresponding to ( |1.1|) p. In this particular case, lamellae 



parallel to the cooling front are first formed and then break up into a triangular 
pattern. In the steady state, a set of three modes, Wi = W^i,25 W^i.s}, 

where each is parallel to one of the three edges of the triangle, invades the 
disordered region. In this illustration, the mode parallel to the cooling front, 
Wi^i, invades first, followed by the remaining two. Under the same bound- 
ary condition, but with different polymer parameters, sometimes a triangular 
lattice is formed by the invasion of the set VV2 = {W^2,i, W^2,2, ^^2,3}, which is 
rotated by 30 degrees with respect to Wi. In general, prior to the establish- 
ment of a steady state 3-mode invasion, there is a competition between Wi 
and W2 (and any other modes which happen to be present). The time evolu- 
tion of the invasion is governed by a set of simultaneous semilinear parabolic 
equations of the formjQ 

dty:}^ = DiA^, + Fi{^i,---,yDN), (i = 1, ■ ■ ■ , iV), (1.2) 

where (pi denotes the amplitude of the i-th mode, is the total number of 
relevant modes, Di is the diffusion constant for the i-th mode, and Fi is the 
'reaction term' (a smooth function). 

In this paper, we first wish to discuss the front selection problem for (|1.2|) : 
when many stable propagating fronts are allowed by the model, what front can 
we actually observe under an ordinary experimental setting? 
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This question is, however, only the starting point of the present pa- 
per, whose main aim is to discuss and illustrate the fundamental role of 
renormalization-group ideas in macroscopic physics. 



The above question about selection has led us to the structural stability 
analysis |ll2[| of ( |1.2|) (Section III). A renormalized perturbation approach is 



given as an algorithm to check the observability criterion due to the structural 
stability analysis (Section IV).[|13| This analysis leads us to a vast frontier of 
renormalization group theory (Section V).|0 In Section II we give a brief 
review of the selection problem. The last section contains a summary and 
comments. This article contains some pedagogical material to clearly demon- 
strate our points, but its main purpose is to announce an intimate relation 
among structural stability, renormalization and singular perturbation. More 
accurate and detailed statements will be published elsewhere. 

II Selection Problem 



The simplest case of ( |1.2| ) is obviously the following scalar equation 

d,^ = dlcp + Ficp). (2.1) 

Fisher introduced the equation with F{{p) = {p{l — ip) (Fisher's equation). We 
assume F(0) = F(l) = 0. If we also assume that Fi^ip) > Vy^ G (0,1), 
then there exists a stable traveling wave solution interpolating between 1 and 
with propagating speed c for all c G [c*,+oo). If F is differentiable at 0, 



then c* > c = 2JF'{0). Thus there are uncountably many stable propagating 



wave solutions for ( |2.1j ). However, usually only one of these is reproducibly 
observable in actual or computer experiments. Thus we have the selection 
problem: what stable traveling wave solution of (|2.1j ) is actually observed? 

To study the selection problem, we must carefully distinguish between 
the model and the system being modeled. We use the word 'system' to denote 
an actual physical system on which we can perform actual experiments. In 
contrast, a model is a mathematical procedure (or equation) describing the 
behavior of some observable(s) of the system which the model is to simulate. 
For example, the model ( |2.1|) simulates front propagation phenomenon such as 
the spreading of an allele of a gene locus in a population (the system). While 
the system apparently exhibits reproducibly a unique propagating front, the 
model allows uncountably many such fronts to exist. What is the selection rule 
for the propagating front solution which corresponds to the actually observed 
front in the system? This is the precise statement of the selection problem. 

In an actual front propagation experiment, say, fire propagation along a 
fuse, we must prepare an initial condition. Fire is set by elevating the fuse 
temperature in front of the observer/experimenter. Thus, in practice the initial 
condition for the system is modified only on a finite region of the system. In the 
model, we must prepare the corresponding initial condition to have a compact 
support. Let us call such an initial condition a physical initial condition. We 
define the 'physical observability' (in the present context) of a solution to 
a given model as follows. If the traveling wave solution is attainable as an 
asymptotic state of the initial value problem with a physical initial condition, 
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we call the traveling wave solution physically observable. This is sensible, since 
we cannot manipulate infinite space to prepare an initial condition. We can 
only modify the system just in front of us. 

Aronson and Weinberger proved the following: 
Theorem A [Aronson and Weinberger |T5|]. For (|2ll ) if F(0) = F(l) = 0, 
F{x) > for any x G (0, 1), and if F'{0) > (these conditions will henceforth 
be called the AW condition) , then the boundaries of any level set for the value 
in (0,1) of the solution with a physical initial condition asymptotically travel 
with the speed c*. □ 

This implies that under the AW condition, the propagating speed we can actu- 
ally observe is the minimum stable speed. We may call this the minimum speed 
principle. Empirically, this is what seems to be generally believed. Certainly, 
we do not have any counterexample for ( p.l| ), even without the AW condition. 
We do not, however, know any rigorous result other than this theorem. 



There is a hypothesis of marginal stability due to Langer. [[161 The linear 
marginal stability analysis is motivated by the following observation. Suppose 
a small localized perturbation is added to the if = state. Since this state 
is unstable, the disturbance grows, and consequently its fronts propagate in 
both directions. We wish to observe the front from a moving frame. If the 
speed of the frame is too slow, the disturbance front outruns us, so that we 
observe a growing disturbance and conclude that = is unstable. If the 
speed of the frame is too fast, we outrun the disturbance, and we say ip = 
is stable. However, the natural front should be self-sustained; the growth of 
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the invading disturbance into the unstable state should be the cause of front 
propagation. Hence, the speed of the front should be the one which makes the 
if = state marginally stable. 

In the moving frame with speed c, ( ^.1| ) reads 

dt^ = {dl + cdi:)ip + F{ip), (2.2) 

where ^ = x — ct. We study the stability of the tip of the traveling wave in 
the following form 

^ = e{t)e'^, (2.3) 
where e is assumed to be very small. We get 

e'{t)=a{k)e{t), (2.4) 

with 

a{k) = k'^ + ck + F'{0). (2.5) 

The marginality condition is Rea{k) = and da{k)/dk = 0. From these, we 
conclude that c = 2^F'{0) is the selected speed according to the hypothesis. 
Notice that this value is the lower bound for the minimum stable speed c 
allowed to the model. Mathematically, we classify ( p.l| ) into two cases|T^: If 
c* = c, the model is called a pulled case, and if c* > c, a pushed case. The 
linear marginal stability analysis works only when the model is pulled. There 
is no established method to distinguish pulled cases from pushed cases. 
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Ill Structural Stability 



To motivate our approach to the selection problem, we first wish to reflect 
upon what we should mean by a good model of a natural phenomenon (or a 
given system). 

Suppose we repeat the same experiment many times and collect data on 
the same observable for a given system. If the observed data cluster around 
some definite value, and the fluctuation around this value is small, we may 
say that the observable is reproducibly observable. Fluctuations around its 
most probable value are due to factors we cannot control. For example, they 
may be due to details in the initial condition or in the system preparation or 
maintenance itself. Now, let us assume that we have a mathematical model 
M of the system under study. If this is a good model of the system, then its 
behavior (at least that corresponding to the reproducible observables) must 
be stable against its modification. That is, in a certain sense, M is close to 
M + 5M, where 5M corresponds to the details beyond our control. 

This is exactly the idea of 'structural stability' of a model first introduced 
in the context of dynamical systems by Andronov and Pontrjagin.[0 Since 
the coefficients of most differential equations important in practice (in physics, 
biology, engineering, etc.) cannot be determined exactly, it is crucial that 
their global features be largely unaffected by tiny changes in these coefficients. 
Therefore, Andronov and Pontrjagin proposed that only structurally stable 
models are good models to do scientific work. An epoch making theorem was 
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later proven by Peixoto |T9| : The set of all the structurally stable C^-vector 



fields on a C°° compact 2-manifold is open and dense in the totality of C^- 
vector fields. This was a very encouraging result, suggesting that we might 
dismiss all the structurally unstable models from science, as suggested by the 
original proposers of the concept. I^C] However, soon it was recognized that if 



the dimension of the manifold is larger than 2, the structurally stable vector 
fields are not dense. 



What does this mean to science? It means at least: 

(i) The World is full of systems which are in a certain respect unstable and 
whose observable results are at least in part irreproducible. 

Then, probably 

(ii) The conventional definition of structural stability is too restrictive for 
science, since the fact that many things are not reproducible is reproducible. 

If there are unstable or irreproducible aspects in the actual system being 
modeled, then a good mathematical model of the system must have features 
unstable with respect to the perturbation corresponding to that causing insta- 
bility in the actual system. Thus a good model should be structurally stable 
with respect to the reproducibly observable aspects, but must be unstable with 
respect to the hard-to-reproduce aspects of the actual system. 

Let us consider Fisher's equation 

dt^ = dl^ + ^{l-^). (3.1) 

We wish to add 6F to its 'reaction' term. If SF is C^-small, that is, \SF\ 
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is small and \SF'\ is also small in [0, 1], then c* changes only a little, and it 
is easy to demonstrate that actually all aspects of the model are structurally 
stable. That is, all changes are continuous with respect to the C^-norm of 6F. 
Unfortunately, it is easy to demonstrate that (|3.1|) is not stable against certain 
C'-perturbations (i.e., without the smallness condition of Consider a 

small spine-like perturbation near the origin. Its size can be made indefinitely 
small while simultaneously making the slope of 6F indefinitely large. Hence, 
we can indefinitely increase the slope of the reaction term at the origin with 
indefinitely C°-small perturbations. This implies that the lower bound c of 
c* can be increased without bound. Hence, the model cannot be structurally 
stable. 

Is this an artifact of the mathematical model and thus a mere pathology? 



Consider the following analogy for (2.1). We may regard the equation to 



be describing the propagation of a flame along a fuse. In this analogy, (f is 
the temperature; is the flash point of the fuse and 1 the steady burning 
temperature. The reaction term F may be regarded as the generation rate of 
heat due to burning (actually, it is the net rate of heat deposition on the fuse: 
the heat generation due to burning minus the loss of heat to the environment. 
In the steady state these must be the same, so F{1) = 0). For ^ 0, we may 
linearize ( p.l|) as 

d,^ = dl^ + F'(0)¥.. (3.2) 

If we put a very small amount of explosive powder along the fuse, we can 
increase -F'(O) considerably. The explosive burning near temperature will 
therefore trigger a very fast propagation of fire along the fuse. Thus, we can 
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imagine an actual system in which a drastic change of c* is possible with a 
very small change of F. We may conclude that the structural instability of the 
model is a desirable feature of a good model. This example thus provides 
an illustration of assertion (ii) above. 

To relax the structural stability requirement of Andronov and Pontrja- 
gin, which requires every aspect [0]of the model to be structurally stable, we 
must consider two things. First, we must require the stability of the model 
only against structural perturbations corresponding to perturbations of the 
actual system which affect its reproducible observables only slightly. We call 
such perturbations physically small perturbations of the system and the cor- 
responding mathematical expressions p-small perturbations of the model. We 
require the structural stability of the model only against p-small perturbations. 
Secondly, we need not require every aspect of the model to be stable against 
p-small perturbations; we have only to require the stability of reproducibly 
observable features. 

Our general conjecture is: solutions structurally stable against p-small 
perturbations describe reproducibly observable phenomena. More precisely, we 
conjecture a structural stability hypothesis: For a good model, only structurally 
stable consequences of the model are reproducibly observable. We must admit 
that there is potentially a tautology here. If we could reproducibly observe a 
phenomenon of a system which is not structurally stable in the model, or if we 
could not reproducibly observe something which the model says is structurally 
stable, then we conclude that the model is not a faithful picture of the system. 
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IV Structurally Stable Solutions of Semilinear 
Parabolic Equations 



For semilinear parabolic equations, we say a C"^-small perturbation is p-small 
if 

sup ^</(||5F||o), (4.1) 

?iG(0,l] 

where 1 1 ■ | |o is the C°-norm, and / is a continuous function such that f{x) 
as X — 0. Notice that the condition has no absolute sign, and only the upper 
bound of 6F/u is specified. Thus we are not demanding the differentiability 
of 6F. 

Now, we have the following theorem: 
Theorem B [Paquette and Oono||T2[] For (^) with F(0) = F(l) = 0, let c*{F) 



be the minimum traveling wave speed for the reaction term F. Then, if 6F is 
p-small, \im\\sF\\,^oc*{F + 6F) = c*{F). □ 



An intuitive idea behind Theorem B is as follows. Suppose ip{x, t) = 0(^) 
(where ^ = x — ct) is a traveling wave solution to ( p.l|) . (p obeys 

f.c|..„)^0, ,4.2) 

or replacing with q, 

q = p, 

dV , , 

P = -cp- (4.3) 
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where F{q) = dV/dq. That is, the problem can be interpreted as a particle 
of unit mass (position q and velocity p) sliding down a potential hill V with 
friction constant c. Hence in this particle analogy, the speed in the original 
problem corresponds to the friction constant. 

A propagating front connecting 1 and corresponds in the particle analogy 
to an orbit connecting the saddle S and the sink (at the origin) O, as shown 
in Fig. 2. If c is too small, the particle overshoots and goes into the region 
q < 0. The corresponding solution of the original partial differential equation 
is thus unstable in the ordinary sense of this word. As can be seen from Fig. 
2, c* is the boundary between overdamped and underdamped motion. Now 
let us put a small potential bump at the origin; this can be done with an 
indefinitely C°-small perturbation to F (or indefinitely C^-small perturbation 
to V). Obviously overdamped saddle-sink connection orbits no longer exist. 
That is, all the front solutions with speed faster than c* are destroyed by this 
perturbation. Obviously, sufficiently underdamped orbits still overshoot the 
origin, so that there must be a boundary between over and underdamped orbits 
which is not far away from the original c*. For c < c*, an appropriate bump 
would convert this c into the critical damping factor (that is, the minimum 
speed of the stable stationary front). However, in this case we can always 
choose a much smaller bump to leave c as an insufficient friction constant 
for the particle to stop at the origin. Hence, the boundary between over 
and underdamped cases must be infinitesimally close to the original c*, if the 
perturbation is infinitesimally small. 
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This intuitive demonstration is technically not easy to rigorize, since al- 
lowed perturbations are not necessarily a simple bump. Still, it captures the 
salient physics (and mathematics) behind the structural stability of c*. 



If g = is not an isolated minimum of V, the propagating solution of (Zl) 



is unique. This can easily be seen from the particle analogy above. Notice that 
it is always possible to eliminate the isolated minimum at g = with a p-small 
perturbation. This, together with Theorem B, implies that c* and only c* is 
structurally stable against physically benign perturbations. 

In the present context, we accept that semilinear parabolic equations are 
good models of front invasion into unstable states. Then the structural sta- 
bility hypothesis implies that the physically observable front speed is the min- 
imum stable speed. If the equation satisfies the AW condition, this is true 
thanks to Aronson and Weinberger's Theorem A. But Theorem B is valid 
even without this condition. 



For the multimode case of ( |1.2[ ), if Fi = ai(pi+ higher order terms, that 
is, if the (f are linearly decoupled, then we can prove a theorem analogous to 
Theorem B. In this case, however, structurally stable speeds need not be 
unique. Generally speaking, there is no further principle to select one among 
the structurally stable speeds. We believe that what we can observe in these 
cases depends on the initial condition. That is, only history can select the 
realized front among the structurally stable ones. Such examples have already 
been given, and in fact, the block copolymer model is one of these. 
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We have been unable to prove the general case where no linear decoupling 
assumption holds for (|1.2|) . Still, we believe that what we have seen for the 
decoupled case holds here too. That is, what we can observe are structurally 
stable fronts, and only history can select the actually realized one among these. 

Why does structural stability imply the minimum speed in this case? The 
key observation to explain this is that the speed c > c* is determined by the 
tip, while the speed c* is determined by the bulk of the propagating front. 
The former may not be hard to understand, because to realize a speed faster 
than c*, we need a fine tuning of the decay rate of the initial condition at 



infinity, as has been demonstrated in the pulled case by Kolmogorov et al. |24 



and KametakapSj. For the pushed case, see [^. The assertion that c* is 



determined by the bulk may sound strange in the case of a pulled front, but 
it is easily seen that even in this case, c* is insensitive to the tip. In both 
the pushed and pulled cases, note that if the initial condition is confined to a 
compact set, or decays to zero more quickly than any exponential, the resulting 
solution decays to zero more quickly than any exponential for all time. Also 
note that if the initial condition decays as ~ exp(— A;a;), where k is at least as 
large as k* (here exp(— A;*x) is the asymptotic form of the steady state solution 
with speed c*), then this asymptotic form is maintained for all time. In all 
of these cases, the asymptotic speed is c*. The initial decay rate therefore 
determines the tip shape for all time, and hence this tip shape has nothing to 
do with the selected speed. Hence, the words 'pushed' and 'pulled' may both 
be misleading. (See [|T2[ for a more detailed explanation.) 
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Now it is easy to understand why the minimum speed is structurally 
stable. Since the tip is extremely fragile against small modification of F near 
the origin, all speeds c > c* are unstable structurally. On the other hand, c* is 
determined by the bulk of the propagating front, which is obviously insensitive 
to a small perturbation. In terms of the fuse analogy, imagine we put a thin 
film of water on the fuse. This would be sufficient to kill the fast propagation 
of fire determined by the tip even if such propagation could be realized in the 
unperturbed system. Thus the structural instability of faster solutions is an 
actual phenomenon; that is, it is not an artifact of the modeling process. In 
this sense, the reaction-diffusion equation is a very good model of, e.g., the 
invasion of a stable phase into an unstable phase. 

Since unstable states are unstable against spontaneous fiuctuations due 
to, e.g., thermal fiuctuations, it is not possible to prepare a wide unstable 
phase region. This is why the moving cooling front is used in the diblock 
copolymer example at the beginning of this paper. Therefore one might think 
that the nonuniqucncss of the propagating front in the model is due to an 
excessive idealization of the actual system: the unstable state of the model 
is really a metastable state with a very small 'activation barrier'. One might 
conclude that this is the reason why in the actual system there is only one 
propagation speed which we observe. We need not deny that there are such 
cases, but in many actual examples, the unstable states are really unstable 
against some particular invasion mode, although they are metastable against 
spontaneous fiuctuations. 
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Consider Fisher's original example of the spreading of an allele in a pop- 
ulation. Of course, the invading allele could be produced de novo by mutation 
in the population, but this is extremely improbable, so the initial population 
is quite stable against spontaneous fluctuations. If the allele is advantageous, 
then the initial population is unstable against its invasion. In the case of the 
fuse analogy we have been using, the flash point Tp is the temperature at 
which the fuel becomes unstable against the invasion of radicals, while the 
ignition point T/ is the temperature at which the fuel can spontaneously pro- 
duce radicals (reacting with oxygen). That is, between Tp and T/, the fuse 
is unstable against the invasion of fire, but metastable (almost stable) against 
spontaneous thermal fluctuations. The distinction between flash point and 
ignition point parallels the distinction between the secondary and primary nu- 
cleation processes. For example, a melt below the melting point should not 
be considered a metastable state when a a crystal nucleus is already present. 
The melt is really unstable against the invasion of the crystal phase. Thus, the 
structural stability requirement cannot be regarded as simply an augmenting 
or auxiliary rule to make excessively idealized models realistic. 



V Renormalization and Structural Stability 

Renormalization group (RG) methods are generally interpreted as a means to 
extract structurally stable features of a model |2^; the structurally stable 
features of the model characterize the universality class to which it belongs. In 
RG terminology Theorem B implies that p-small perturbations are marginal 
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perturbations for c*, but that some p-small perturbations are relevant to speeds 
larger than c*. Furthermore, we know that generally speaking, C°-small per- 
turbations could be relevant. 

Thus Theorem B affords a method to judge whether the front with speed 
Co is observable or not through the study of its response to 5F corresponding 
to a small potential bump added to the model: if the change of the speed 5c 
vanishes in the limit of vanishing bump (that is, if 5F is a marginal perturba- 
tion), then Co is observable. Otherwise, cq is not observable. As we have found, 
this procedure works numerically. In response to a p-small perturbation 5F, 
the change in the speed of ( |1.2| ) observable in numerical computations vanishes 
with ||(5F||o. Let us consider an example. 

As stated above, we have been unable to prove a statement analogous 
to Theorem B for multi-mode equations which display linear order coupling. 
We believe, of course, that our structural stability hypothesis applies to these 
equations as well, and in support of this conjecture, consider one such model 
for the present study. We note that similar behavior can also be easily observed 
for single-mode and multi-mode, linearly decoupled equations. Consider the 
following model equation: 

d^i/dt = 2A^i + Fi(V^i,V^2) 

dij2/dt = AV'2 + F2(^l,V'2), (5.1) 

where Fi = ipi + \ip2 — ipf, and F2 = ?>ip2 + — i^f - — i^2- We numerically 
studied the behavior of ( p.l|) in response to the perturbation Fi Fi + 6F1 
and F2 ^ F2 + 6F2, where 6Fi = —lOipi if ipi < e and otherwise. Note that 
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{SFi,SF2) can be considered as the discretization of a p-small perturbation. 
{6Fi,6F2) is analogous to the film of water discussed in the context of the 



fuse analogy. If a traveling wave solution of (5J.) with speed c is observable 
(structurally stable), the speed of the observable solution of the perturbed 
equation must converge to c as e — 0. We numerically determined the ob- 
servable propagation speed of the unperturbed equation, as well as those of 
perturbed equations with several values of e. The results of this study, shown 
in Table I, support our structural stability hypothesis; the observable speed 
changes continuously in response to a p-small perturbation. 

We next studied a "tip driven" solution (as opposed to the "bulk driven" 
solution considered above) of (|5.1|) . We were able to produce such a solution by 
choosing two small positive values ei and e2, and forcing the value of x at which 
both ipi = ei and ip2 = ^2 to move at speed c = 10. We chose ei = 0.248 x 10^^^ 
and €2 = 10~^^. With these values, the eigenfunction of the linear equation 
corresponding to ( ^.1|) for the traveling wave solution with c = 10 is given by: 
const. [ei, 62] exp(—A;x), with k = 0.323. We then computed the speed of the 
resulting front by watching the point at which ipi = 0.01. Not surprisingly, 
this value was 10. However, when we applied perturbations to the tip driven 
model identical to those applied to the bulk driven model, in each case, the 
propagation speed computed was also identical to that found for the bulk 
driven model. For the tip driven solution, the response of the model remains 
finite as the size of the perturbation vanishes. The above considerations thus 
lead us to conclude correctly that it is unobservable. 
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Once more returning to the propagation of fire as a physical analogy, this 
result can be interpreted as follows. For the dry fuse, we are able to force 
the system to exhibit 'fast' flame propagation by running a torch along the 
fuse to ignite it at the desired speed. When we add a film of water which the 
torch is not able to evaporate as it runs past, however, the behavior of this 
torched system cannot be distinguished from that of the untorched system. 
Its response to this small perturbation is therefore large. 

Let 00 be a stable traveling wave solution of (|2.1| ) with speed cq. Let us 
add a p-small structural perturbation 5F to ( p.l|) with ||5-F||o of order e, and 
assume that in response the front solution is modified to 0o + ^0- Linearizing 
( p.l|) to order e in the moving frame with velocity Cq, we obtain formally the 
following naive perturbation result: 

rt r+oo 

<50(e, t) = e-^°«/2 / dt' / t; e', t')e'^'^ ^'SFiMO)- (5-2) 

J to J ~oo 

Here to is a certain time before 5F(0o(O) becomes nonzero, and G is the 
Green's function satisfying 

f)C 

__CG = 6{t-m^-a (5.3) 
with G — > in 1^ — oo, where 

^-^ + ^00(0)- 1- (5.4) 

Since by C°-infinitesimally modifying F, we can always cause C to have 
as an isolated eigenvalue, we may safely disregard all possible complications 
introduced by the presence of a eigenvalue which is not isolated from the 
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essential spectrum. Formally, G reads 

G{^, t- e', t') = uoiCKiO + E e-^"(*-*')«„(0<(e'), (5.5) 

where Cuo = 0, and Cun = XnUn- The summation symbol, which may imply 
appropriate integration, is over the spectrum other than the point spectrum 
{0}. Since the model is translationally symmetric, Uq oc e'^''^/^0Q(^). Due to 
the known stability of the propagating wavefront, the operator C is dissipative, 
so is the least upper bound of its spectrum. Hence, only uq contributes to 
the secular term (the term proportional to t — to) in ^'P- Thus we can write 



6(t>=-{t-to)5c4>'o{0 + mn (5.6) 



where the suffix B means "bare', {5(j))r is the bounded piece (regular part), 
and 

One may immediately guess that this 5c is the change in the front speed, but 
the naive perturbation theory is not controlled. A renormalization procedure 
can be used to justify the guess as follows. [O 



The first term in (|5.6|) is divergent in the ^ -co limit. We introduce an 
arbitrary subtraction factor /i to separate the divergence by splitting t — to as 
t — fi — (to — fi) , and then absorb the divergence — to through renormalization 
of 0o(O to 0-r(^7/^). To order e we get 

(Pb{0 = 0fl(^> /^) - ^c{t - /^)0i?(^, /^), (5.8) 

where 0o in the second term is replaced with 0^, because 6c is already of order 
e, as seen from (|5.7|). The RG equation is 90B(^)/9/i = 0. Hence, to order e 
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the RG equation is, after equating /i with t, 



^ + 6c^ = 0. (5.9) 



dt 

Thus the speed of the renormahzed wave is indeed Cq + Sc. 



The formal expression (|5.7| ) is legitimate only when both 5F and —5F are 
p-small. That is, the formula is legitimate only when 6F is linearizable near 
the origin. Since we do not know whether the renormahzed perturbation result 
is asymptotic or not, strictly speaking the formal expression ( |5.7|) and the true 
change 6c = c{F + 6F) — c{F) itself should be distinguished. Furthermore, 
the expansion is correct only if the terms obtained are finite, so if c is not 
structurally stable, the formal expression may not be justified. Still, (|5.71 ) 
seems to give us the correct information about the observability of c. 

For example, if we add 6F = e0(l — (p) to ( |3.1D with F = 0(1 — 0), then 
( ^.7|) gives c* ^ 2 + e; the exact result is, of course, c* = 2^1 + e. If we add 
SF = 9{(f) - A)(0 - A)(l - 0) - 0(1 - 0), with A > and 6* being the unit 



step function, then 6c = vA if cq = 2, and 6c = ycg — 4 in the A — > limit. 
Hence, only when Cq = 2 does c change continuously with the perturbation. 



VI Singular Perturbation and Renormalization 

The reader may make the criticism that the renormalization approach in the 
preceding section is nothing but a singular perturbation approach (the method 
of stretched coordinate). Why do we need such a (purportedly) heavy ma- 
chinery as RG? Before answering this question, we must stress that RG is not 
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an esoteric machinery. As mentioned in the preceding section, it is a (the?) 
method to extract structurally stable features of a given model. For example, 
in the case of critical phenomena, we wish to study global features which are 
insensitive to small scale details. That is, we are pursuing the features of the 
model stable against structural perturbations corresponding to the small scale 
details. 

In this section, we first demonstrate that the calculation in the preceding 
section is just the standard renormalization group theory for partial differential 



equations. [p!4| , P8| Then, we demonstrate that the ordinary singular perturba- 
tion method is understood very naturally from the RG point of view. Actually, 
we wish to claim that singular perturbations are most naturally understood 
as renormalized perturbations. 

Introducing new variables X = and T = e*, the propagating front so- 
lution reads (f){x — ct) = $(XT~^). Thus the front speed is interpreted as an 
anomalous dimension. This is obvious; since the variables inside logarithms 
must be dimensionless, c cannot be determined by dimensional analysis. If 
we introduce Tq, defined by to = IhTq, then t — to = ln(T/To). From this we 
may interpret Tq as an "ultraviolet cutoff" scale. Hence, the to ^ "Oo limit 
corresponds to the cutoff limit in the usual field theoretic calculation or 
in our PDE calculation. In the ordinary multiplicative renormalization group 



scheme, |29] the logarithmic singularity ln(T/To) is absorbed into the renor- 
malization group constants. Usually, we introduce an arbitrary length scale 
L and rewrite T/Tq as (T / L){L/Tq). ln(L/To) is then removed by renormal- 
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ization. Our /i above is nothing but InL, and the sphtting of the logarithmic 
terms should correspond to the splitting t — + — to- ^ represents the 
divergence to be absorbed into some phenomenological parameter. 

Now, with the aid of the presumably simplest (but representative) exam- 
ple, we demonstrate our point that singular perturbation is best understood 
as renormalized perturbation. Consider the following linear ODE: 

ex + X + X = 0. (6.1) 

We pretend that we cannot obtain its closed analytic solution and apply a very 
simple-minded perturbation approach. Expand x as x = xq + exi + ■ ■ ■. We 
have 

xo + xo = 0, (6.2) 
xi + xi = —xq. (6.3) 

Solving these equations, we can easily get the following formal expansion: 

X = Aoe-(*-*o) - eAo{t - to)e-(*-*'') + 0[e''], (6.4) 

where Aq is a constant determined by some initial condition. Now, the second 
term contains the prefactor t — to, and is thus a secular term; the ratio of 
the first and the zeroth order terms diverges in the to ^ —oo limit. As done 
above, we now introduce fi, split t — to as t — fi + fi — to, and absorb fi — to 
into Ao, which is due to the initial condition we do not know. In this way, Ao 
is renormalized to A. We rewrite (|6.4| ) as 

X = Ae-^'-^^ - eA{t - /i)e~(*-^) + 0[e^]. (6.5) 
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Since /i is not in the original problem, obviously dx/d^ = 0. This is the 
renormalization group equation. After differentiating ( |6.5|) with fi and then 
setting fi equal to t, we get 

^ + A + eA = 0[e% (6.6) 

This is exactly the equation obtainable, for example, by the reconstitution 
method pO|. Solving this equation (ignoring the second order term), and 
putting the result into ( |6.5|) with /i = t, we get 

X = 5e'(^+^)*, (6.7) 

where B is the 'phenomenological constant' we must fix appropriately to re- 
produce the observable result. Clearly ( |6.7| ) is the formula obtained by the 
usual stretched coordinate method, or a multiscale expansion scheme. Here 
the result is obtained without the introduction of modified variables or coor- 
dinates. 

One might think this agreement is only fortuitous. To see that this is not 
the case, consider ( |6.4| ) again. This formula is reliable if e(t — to) is sufficiently 
small. Instead of calculating the result at t at once from to, we could proceed 
step by step just as in the Wilson renormalization group theory. Let us 
divide t into time spans and first solve the problem from to t/N (for 
simplicity, we set to fo be 0). We get 

x{t/N) = 5e~*/^(l - et/N) + 0[{e/Ny]. (6.8) 

Now use this as the initial value and solve x{2t/N) to order e, etc. We even- 
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x{t) = B [e-*/^(l - et/N)] ^ . (6.9) 



tually get 



Taking the N —>■ oo limit, we get 

To obtain a solution reliable not only for large t but for all t, in the 
standard singular perturbation procedure, the so-called inner and outer ex- 
pansion and their matching are required. [^] Now we demonstrate that from 
only the inner expansion, we can construct a uniformly valid solution by a 
renormalization group method. 



First is rewritten as 

x" + x' + ex = 0, (6.10) 

where ' implies the derivative with respect to r = t/e. Naive perturbation 
gives the following result: 

x = Ao + Boe~^ - e[Ao{T - 1 + e^O + 5o(l - re'^ - e~")] + 0[e^]. (6.11) 

Introducing fi into the secular terms through t ^ t — fi + fi, we wish to absorb 
H (here tq is set to be by an appropriate time shift) by renormalizing Aq 
and Bq. Let us proceed more systematically by introducing the multiplicative 
renormalization factors, = 1 + eai + ■ ■ ■ and Zb = 1 + ebi + ■ ■ ■, and 
renormalized coefficients as A = ZaAq and B = ZbBq. Putting everything 



into ( |6.11| ), we get to order e 

X = A{l-eai-\ ) + B{l-ebi-\ )e"^ 

-e[A(r - fi + + e~^) + 5(1 - (r - /i + fi)e-^ - e"^)] 

+0[e% (6.12) 
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Thus the choice ai = jj, and hi = —jj, successfully ehminates the secular terms, 
and we get the renormalized perturbation result 

X = A + Be-^ -e[A{r-^-l + e-^) + B{l-{T-iJ.)e'^ -e-^)]+0[e% (6.13) 

Notice that A and B are now functions of /i. Since x should not depend 
on /i, which is introduced independent of the original problem, we have the 
renormalization group equation dx/dfi = 0. From (|6.13| ) we get 

= ^ + ^e-^ - e[-A + Be--] + 0[e']. (6.14) 

Here we have used the fact that derivatives are of order e. Due to the functional 
independence of 1 and e~'^, we get 

dA , dB , , 

— = -eA, — = +tB. 6.15 
djjL djj, 

Solving these and equating /z and r in (|6.13|) , we get 

xr = Ae-'- + Be-^^-'^- + e{A - B){1 - e~"). (6.16) 

Let us compare this with the result obtained by the standard inner-outer 
matching method to order e (that is, both the inner and outer solutions are 
obtained to order e; notice that this calculation is partially second order): 

X = Ae-'- + Be-- + Bere^ + e{A - B){e-'- - e"") - eMre"^". (6.17) 

Except for the term, all the terms are correctly given by the RG procedure. 
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VII Reductive Perturbation and Renormaliza- 
tion 



Now, let us look at ( ^.91) . This is the equation of the wavefront as seen from 
the moving coordinate translating with the speed of the unperturbed front. 
From this frame the motion of the perturbed front is very slow. Hence, (|5.9| ) 
is regarded as a slow-motion equation, like an amphtude equation obtained by 
the so-called reductive perturbation methods. That an amplitude equation 



is an RG equation is not a fortuitous relation but a rule. 

To see the point, let us consider the following slightly dissipative nonlinear 
hyperbolic equation: 

where A(m) is a sufficiently smooth function of u, and 77 is a positive constant. 
We consider a small amplitude wave in the background of the constant solution 

Mo, 

u = uq + eui + e^U2 H , (7.2) 

where e denotes the amplitude of the wave. 

First we study the case without dissipation (77 = 0). Let us perform a 
naive perturbation approach. Let Aq = A(mo). We have 

dtUi + XqO^Ui = 0, (7.3) 
dtU2 + \od^U2 = -\'{uo)uid^ui, (7.4) 

and so forth. Introducing independent variables (^ = x — X^t, t) to replace 
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(x, t), these equations can be rewritten as (notice that dt now reads dt + Xodx) 

dtux = 0, (7.5) 
dtU2 = -\\uo)uid^ui. (7.6) 

Thus the right hand side of the second equation is a function solely of ^ in this 
coordinate system, so that it gives a secular term. Thus to order we have 
the following general solution: 

u = uo + eFo(0 - e\t - to)A'(Mo)i^o(0^o(0- (7-7) 

We introduce /i as we did for the propagation wave and split t — to as t — yU + 
/i — to- Then we absorb — to into the renormalized version F(^, jj) of -fo(0- 
The renormalized perturbation result reads to order 

u = uo + eF(e, /i) -e\t- fi)X'{uo)F{^, f^)d^F{^, fi). (7.8) 

The renormalization group equation must be du/d^ = 0, so that we get to 
order 

df,F + e\\uo)Fd^F = 0. (7.9) 
If we identify and t in (ItTsI) , we get u = uo + F{$,, t), so ( [7. 91) with /i = t, or 

+ Aoa,F + eX'iuo)Fd.,F = 0, (7.10) 

in the original coordinate system is the equation of motion for the small am- 
plitude wave. 

With the introduction of a weak dissipation, the first equation of (|7.6| ) 
should not be affected (this is the precise meaning of weak dissipation). At 
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worst, only the second equation is modified as 



dtU2 = -X'{uo)uid^ui + {r]/e)d'^ui. 



(7.11) 



Thus (fTTp is modified to be 



u = uo + eFo(0 -it- to)[e'A'(Mo)Fo(OF^(0 + evdiFoiO]- (7.12) 



This is of course a standard resuh obtained by a reductive perturbation method. 



At the beginning of this paper, we illustrated how to generalize the concept 
of structural stability so that it is not excessively restrictive, and we applied 
it to the selection problem of front propagation speeds. Since the basic idea 
of renormalization group theory is to extract structurally stable features of 
a given model, this consideration naturally led us to the RG study of front 
propagation. 

This study in turn revealed two very general conclusions, which are illus- 
trated with simple examples: 

(1) Singular perturbation methods are best understood as renormalized per- 
turbation methods, and 

(2) Amplitude equations are just RG equations. 



Hence instead of (|7.9|) , we arrive at Burgers equation: 



hF + e\'{uo)Fd^F - 7]dlF = 0. 



(7.13) 



VIII Summary 
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The latter in particular strengthens our belief that RG is a prerequisite to do 
physics without being affected by unknown (high-energy) details of the world. 
A more systematic presentation with numerous examples of (1) and (2) as well 
as the relations to the solvability condition, center manifold theory |Q, etc. 
will be given elsewhere. 
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Table Caption 

Table I. The observed speed of the front as a function of the size of per- 
turbation. The speed is a continuous function of perturbation. That is, this 
observable speed is structurally stable. 

Figure Captions 

Figure 1. A cell dynamics model simulation (for details, see 0) of a block 
copolymer film with an invading triangular phase. Initially, a periodic pattern 
is imposed on the right edge of the system. As time proceeds, clearly the 
lamellar mode parallel to the invasion front leads the ordering process into the 
unstable uniform phase. Eventually the lamellar pattern breaks up into the 
final triangular pattern (defects may be introduced in this example because 
of a slight mismatching of the parameters and the systm size). Thus, Wi^i 
invades first, followed by the remaining two modes. 

Figure 2. An intuitive explanation of the structural stability of the slowest 
stable propagation speed c*. The trajectories corresponding to the traveling 
wave solutions are illustrated for ( [4.3| ). The left column with U is for the un- 
perturbed model, and the right column with P for the model perturbed with 
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a small potential bump at the origin. S is the saddle, and A is the newly 
formed stable point with the potential bump. The friction constant c (that 
is, the front propagation speed in the original problem) is decreased from A 
to C of the figure for both columns. BU illustrates the critical speed c* case; 
if c is slightly decreased further, then the trajectory overshoots the origin as 
CU. The potential bump at the origin prevents all the overdamped trajectories 
like AU from reaching the origin, as illustrated in AP. For most underdamped 
cases like CU, a small bump is not enough to stop overshooting. Between AP 
and CP there must be a critical friction coefficient for the perturbed model, 
but it must not be far away from the unperturbed one. Hence, c* must be 
structurally stable. Furthermore, no other c can be structurally stable. 
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